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Particle correlation measurements associated with a hard or semi-hard trigger in heavy-ion colli- 
sions may reflect Mach cone Shockwaves excited in the bulk medium by partonic energy loss. This is 
of great interest because, when compared with theory, such measurements can provide information 
on the transport properties of the medium. Specifically, the formation of Mach cone Shockwaves is 
sensitive to the viscosity and speed of sound, as well as the detailed nature of the jet medium inter- 
action. However, modeling the physics of Shockwave excitation to obtain a meaningful comparison 
with the measured correlations is very challenging since the correlations arise from an interplay of 
perturbative as well as non-perturbative phenomena at different momentum scales. In this work 
we take a step in that direction by presenting a systematic study of the dependence of azimuthal 
particle correlations on the spatio-temporal structure of energy deposition into the medium. Our 
results indicate that detailed modeling of the evolution of an initially produced hard parton and the 
interaction of this evolving state with the medium is crucial, as both magnitude and shape of the 
Shockwave signal show a strong dependence on the assumptions being made. 

PACS numbers: 12.38.Mh 



I. INTRODUCTION 

It is useful to describe the dynamics in ultrarelativistic 
heavy-ion collisions (URHIC) in terms of the 'bulk' and 
of 'probes'. The bulk medium, that is, the quark gluon 
plasma (QGP), describes QCD matter produced in the 
collision which is strongly coupled, exhibits collectivity 
and behaves like a thcrmalizcd, near-perfect liquid. The 
bulk matter is chiefly responsible for the properties of low 
transverse momentum (Pt) hadron production. On the 
other hand, there arc also hadrons produced at larger mo- 
menta which arc clearly not thcrmalizcd. Such hadrons 
originate from hard partonic processes which probe such 
small time and distance scales that they are essentially 
unmodified by the medium. However, these high px par- 
tons undergo a final state interaction while they propa- 
gate through the medium before hadronization. This 'jet 
quenching' [H-Q has been expected by theory and is ex- 
perimentally confirmed in measurements of the nuclear 
suppression factor Raa 0- 

If there is energy loss from a high px parton, energy 
conservation requires that this lost energy flows some- 
where. Measurements of particles associated with a high 
Px trigger hadron [H, Q have given a hint to answer this 
question: Instead of a back-to-back jet-like correlation 
structure on the near (trigger) side and the away side 
as observed in d-Au collisions, the observed correlation 
shows a surprising splitting of the away side peak into 
a double-hump structure. This has been early on inter- 
preted as the reaction of the bulk medium to the hard 
probe in terms of a Shockwave [To| and early phenomeno- 
logical investigations with a fluid-dynamics inspired ap- 
proach to propagate the Shockwave in the background 



of an evolving bulk medium and with a full modeling 
of the trigger bias have confirmed that a shock- 

wave signature is not erased by the medium-flow induced 
distortion or by the averaging over many different trig- 
gered events, but can indeed account for the correlation 
observed in the data. Since then, the focus has been 
on a more rigorous theoretical formulation of coupling a 
source of energy and momentum into the hydrodynam- 
ical equations using Hard Thermal Loop or AdS/CFT 
methods [Hfil. 

In bringing such proof-of-concept calculations closer 
to a comparison with data, one of the key questions is 
the time dependence of the energy and momentum de- 
position into the medium. In this paper, we make the 
assumption that the energy and momentum absorbed by 
the medium are related by an on-shell condition so that 
udE/dt = dp/dt (see Section [II Cj) , where dE/dt and 
dp/dt are the energy and momentum deposition rates, 
respectively, and u is the velocity of the high px par- 
ton. For this reason, we focus on the time dependence 
of the energy deposition, which is potentially driven by 
many different effects: First, the strength of the inter- 
action between hard parton and medium depends on the 
medium density and in a real heavy ion collision this den- 
sity varies as a function of space roughly as given by the 
nuclear overlap and, due to the expansion of the medium, 
also drops as a function of time. The expansion dynam- 
ics therefore tends to lead to less interaction with the 
medium and hence less energy loss at late times. On the 
other hand, radiative energy loss in a constant medium of 
length L has a characteristic L 2 dependence due to LPM 
suppression of near collinear gluon radiation, and this ef- 
fect tends to increase energy loss at late times, to some 
degree even in an expanding medium. In addition, in (l9j 
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it was suggested that gluons radiated from a hard parton 
subsequently themselves interact with the medium and 
hence contribute to the energy flow into the medium, 
leading to a 'crescendo' in the Shockwave excitation and 
large energy deposition at late times. On the other hand, 
low energy partons undergoing strong energy loss cannot 
act for a long time as sources of energy, but become ab- 
sorbed by the medium after a short time already. This 
effect again tends to lead to small energy deposition at 
late times, as energy deposited into the medium early on 
is not available later. 

The arguments given so far assume that the source 
of energy and momentum entering the medium is a sin- 
gle on-shell parton which subsequently undergoes inter- 
actions which lead to induced gluon radiation. However, 
in a typical hard event, partons are produced with large 
initial virtualities and even in vacuum evolve into a par- 
ton shower where the individual quanta have lower virtu- 
alities, and the timescale of the shower evolution is such 
that it at least partially takes place before a medium is 
produced. This implies that it may be wrong to think of a 
single parton initially depositing energy into the medium 
— the dynamics may rather be that a developed parton 
shower acts as a strong source initially, but energy depo- 
sition decreases soon as the energy of subleading shower 
partons is quickly depleted. 

It follows from the above that modeling the time- 
dependence of energy deposition correctly is not a simple 
and straightforward issue. It is the purpose of this paper 
to demonstrate that the question is nevertheless highly 
relevant: Modeling the time-dependence of energy depo- 
sition in a different way alters both the magnitude and 
the shape of the Shockwave signal in a significant way. 
We illustrate this point using different model assump- 
tions for the energy deposition within a constant medium 
and solving linearized hydrodynamical equations. The 
paper is organized as follows: In section III Al we review 
the underlying formalism of linearized hydrodynamics, 
while in section III Bl we discuss how to obtain the az- 
imuthal hadron spectrum from the linearized equations 
of motion. In section III CI we discuss the motivation for 
the form of the hydrodynamic source term used in our 
results and show how it depends sensitively on the time- 
dependence of the energy deposition. In section III Dl we 
present the different energy deposition scenarios used in 
our results and discuss the underlying physics assump- 
tions of each one. The reader only interested in the re- 
sults of our calculations can skip directly to sections [HI Al 
and IIIIBI where specific numerical inputs and resulting 
azimuthal spectrums are presented. 



daries) propagating in an infinite and static QGP. Ad- 
ditionally, we ignore any net baryon density, which is a 
reasonable assumption for RHIC energies at mid-rapidity 
[2l| . This parton acts as a source of energy and mo- 
mentum which is coupled to the linearized hydrodynamic 
equations of the underlying medium. The linearized ap- 
proximation is valid when the energy and momentum 
density generated by the hard parton is small compared 
to the equilibrium energy density of the medium. More 
will be said on the linearized approximation in Section 

EH 

The assumption of an infinite and static QGP is clearly 
unrealistic for heavy-ion collisions. However, our purpose 
here is not to present a study which is directly compa- 
rable to experimental data, but rather to show the ef- 
fect of the time-dependence of energy deposition on the 
azimuthal particle spectrum associated with a hard par- 
ton. The linearized approximation in a static medium is 
a good toy model for such a study because the effect of 
changing parameters such as viscosity or energy deposi- 
tion scenarios can be easily extracted. Even within the 
linearized approximation one could go beyond the static 
medium, for instance by assuming an underlying Bjorkcn 
expansion. This, however, would introduce complications 
associated with underlying flow fields and boundary con- 
ditions, and is beyond the scope of the current paper. 

In the linearized approximation, the effect of the source 
is to create a local perturbation in the medium, so that 
the energy-momentum tensor has the linearized form 



(1) 



where <5T Ml/ is the perturbation generated by the source, 
and Tq V is the equilibrium energy-momentum tensor of 
the underlying medium. The fast parton's ability to per- 
turb the medium is encoded in the source term, J v (to 
be specified below), which couples to the gradient of the 
energy-momentum tensor 



(2) 



where d^Tg" = 0. 

The equations of motion for a medium coupled to a 
source in linearized hydrod yna mics are discussed in sev- 
eral places (for instance, jiol . [l6j). The solution for 
8T^ U in terms of J v is most easily expressed in momen- 
tum space by taking the Fourier transform of The 
result to first order in shear viscosity, 77, for the per- 
turbed energy density, ST aa = 5t, and momentum den- 
sity, 8T° l = g, are given by 



II. THE MEDIUM EXCITATION 

A. Linearized Hydrodynamics 

In what follows we consider a hard parton (parton here 
refers to the parent parton and the associated secon- 



ded, w) 
gi(k,w) 
gr(k,w) 



ik Ji(k, to) + J°(k, u))[iu - T s k 2 ) 
uj 2 — c 2 s k 2 + iT s u>k 2 
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In the above result, c s denotes the speed of sound, 
T s = 57 4 '' — r = 4^ is the sound attenuation length, e 

s 3(e +Po) 3sT o ' u 

and p are the unperturbed energy density and pressure, 
respectively and s is the entropy density. Also, the source 
and perturbed momentum density vectors are divided 
into transverse and longitudinal parts: g = kg^ + g^ 
and J = kJi + Jy, with k denoting the unit vector in 
the direction of k. The position space result for equa- 
tions (0 - [5]) are obtained by reverse Fourier transform 
using the general rule 



F(x,t) 



1 



(2tt)' 



d k / duj e 



ik-x— icut 



F(k,w). (6) 



B. The Azimuthal Spectrum 

Once the source term has been specified, and a solu- 
tion for 5e and g is obtained from equations ([3] -[6]), we 
will be interested in determining the azimuthal particle 
spectrum generated by the source. Ignoring viscous cor- 
rections, the induced medium flow velocity is given by 



Su(x) 



g 



eo+Po e (l + c 2 s ) 



(7) 



where we have used c 2 s = dp/de m po/eo for T ^> T c [22J. 
The total medium four-velocity, which is a sum of the 
underlying medium and induced velocities, is given by 



(m + 8u{x)Y = 1 



g 



eo(l 



(8) 



where Su° = in the limit of a static background. An ex- 
pression for ST can similarly be found from dimensional 
considerations. We write the medium energy density as 
e = AT 4 , where A is some constant, from which one has 



leading to 



e + (5e«^T 4 (l + 4 



ST(x) = ^T . 



— ) 
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(9) 



(10) 



Having expressions for the flow velocity and tempera- 
ture, it is now possible to construct the medium's distri- 
bution function, which in the Boltzmann limit is given 
by 



f(x,p) 



exp 



(uq + 5u{x)Yp ll 



To + ST(x) 



(11) 



Here, j3 = 1/T is the inverse temperature. The distribu- 
tion is converted into an azimuthal particle distribution 
by using a Cooper-Frye freeze-out scenario 12311 . Con- 
sistent with the approach discussed in [13, Il8j . the fi- 
nal azimuthal particle spectrum for massless particles at 



mid-rapidity (y 

dN , 

:(y = o) = 



0) is given by 



dyt 



I't 



>T dprPT 

(2tt) 3 



d^p»(f(x,p) 



where S„ is the freeze-out hypersurfacc and 



(pt,0,Pt sm((j)),pT cos(</))). 



/o) 
(12) 

(13) 



The isotropic background contribution, fo = e~^ 0U ° PlJ , 
is subtracted in (TT2l . In what follows, it is understood 
that the source parton propagates along the z axis, which 
also determines the direction of <f> = in (fTB"]) . We will 
consider an isochronous freeze-out, as appropriate for a 
constant medium, in which case cE^ = dV(l, 0). 



C. The Source Term 

Nothing has been said up to this point about the form 
of J v to be used in ([3]-[5]). A common choice is the simple 
form 



J v (x) = —5(x-ut)W 



(14) 



where U v = (l,u), u is the velocity of the source par- 
ton, which is assumed to be at the origin at t = 0, and 
dE/dt is the time-dependent rate of energy loss into the 
medium (in principle this can be less than the total rate 
of energy loss, if energy lost from a leading parton is for 
example redistributed to non-thermalized degrees of free- 
dom) . This expression contains an cikonal assumption as 
it describes a point source propagating with the speed of 
light on a straight line. The advantage of the source term 
in (|14p is that it trivially conserves energy and momen- 
tum at each step in time. One can see this by integrating 
both sides of ([2]) over all space, from which it is found for 
the v = component 



d I A A dE 

— I dxde = — — , 
dt dt 



(15) 



which shows that the energy going into the medium is 
properly accounted for. A similar exercise holds for the 
other components of v as well. 

However, one can go beyond the simple form of (TT4")) 
and still conserve energy and momentum at each step in 
time. Consider adding a total derivative to (fT4"]) 



dE 
~dt 



{U v - Aa i/ )<S(x-ut) 



(16) 



where A is a coefficient with dimension of length. Again 
integrating both sides of © over all space shows that 
energy and momentum are still properly accounted for. 
A acts as a local medium excitation parameter, that is, its 
contribution integrates to zero globally. The replacement 
in (|16[) is motivated by the form of the kinetic theory 
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derived source term for a parton in a perturbativc QGP 
obtained in [l5| , and is similar to the source term derived 
for a quark in a strongly-coupled supersymmctric Yang- 
Mills plasma (25|. The relativistic limit (7 3> 1) of the 
source derived in can be put into the form 



r 



a s C 2 m 2 - ) 
8tt 



(l,u) 



7 



(p 2 +-f 2 z 2 _) 3 /< 



-d v 



1 



2(p 2 + 7 2 z 2 ) 
(17) 



where G2 is the Casimir of the source parton, mo is the 
Debye screening mass in the medium, p = \/ x 2 + y 2 and 
z_ = z — ut for a source parton propagating in the posi- 
tive z direction and 7 = 1/ \[\ — u 2 . 

If one considers the source to be localized, (|17p can be 
put into the mold of (flU]) by integrating the distributions 
over all space, and normalizing to a <5 function. It is help- 
ful to add a damping factor, e~ prriD , to the distributions 
in (|17[) which regulates an infrared divergence that arises 
when integrating over all space. The form of our damp- 
ing factor is motivated by the fact that medium induced 
screening of the hard parton's color fields occurs on the 
Debye scale. The damping factor simulates the Debye 
screening of a Lorentz contracted distribution, hence the 
form e~P m °. 

We find 



(p 2 +7 2 z 2 _) 3 /2 



4tt Go 



m D 



5(x - ut) (18) 



where (2^E p T)^ 1 has been introduced as a short dis- 
tance cutoff (E p is the energy of the hard parton and T 
the medium temperature), and Go is a representation of 
the incomplete Gamma function 



G (z) 



dt- 



One can also show 



1 



2(p 2 + 7 2 z 2 ) 7TOD 



5{x - ut). 



(19) 



(20) 



With these approximations, (|17[) is thus written in the 
form of ([TBI) as 



J" = ^-{U u -\d u )5{x-ut) 



where 



and 



4 7 m D G 



dE _ a s C 2 m 2 D ^ I m D 
~dt ~ 2 ° 1 2/EpT 



(21) 



(22) 



(23) 



Localizing the source term to a 6 function, such as what 
has been done in (f2"Tj) , has the advantage of making the 
linearized hydrodynamics easier to solve, but also has a 
physics justification. Hydrodynamics is a long distance 
theory which assumes local thermal equilibrium, whereas 
the energy deposition occurs on short distance scales and 
is a highly dissipative process. There is a natural separa- 
tion of distance scales between the energy deposition and 
the medium response as described by hydrodynamics (in 
perturbation theory these scales are l/(gT) and l/(g 4 T), 
respectively, where j< 1). Localizing the source term to 
a d function is consistent with this separation of scales. 

As will be seen in the results below, the coefficient A 
is especially important for generating Mach-like signals 
in the final azimuthal particle spectrum. This can be 
traced back to equations (J3] - [5j) . Equation ([5]) is a diffu- 
sion equation and the quantity g^ is diffusive momentum 
density generated by the fast parton. Physically, the dif- 
fusive momentum contribution is a wake which flows in 
the direction of the fast parton's propagation and is not 
a sound wave. Previous studies [171 124| have shown that 
the diffusive momentum tends to fill up any double-peak 
structure in the final spectrum. On the other hand, equa- 
tions (J3]) and (|4]) describe damped sound waves propagat- 
ing at speed c s : it's clear that Se and g£ are the energy 
and momentum density carried by sound generated by 
the fast parton and will be responsible for Mach-like sig- 
nals in the azimuthal spectrum. When writing a source 
of the form (|16|) in momentum space we find 



/HE 
d 4 xe- lk - x+lult — (t) (U v - A0")<J(x- ut) 
dt 



dE 

dt^—(t) e- <fc « t+&rt (U v + iX k v ) 
dt 



(24) 



where a time derivative on dE/dt as well as boundary 
terms have been ignored, and the source is assumed to 
propagate from time t = to t = T. We will discuss why 
we have dropped the derivative and boundary terms at 
the end of the results section III B. Here we simply note 
that the derivative term is numerically insignificant for 
the energy deposition scenarios we consider, and that the 
boundary terms arc an artifact of stopping and starting 
the source at a specific moment of time and obscure the 
physics we are trying to explore. It is immediately clear 
from (pM)) that the term proportional to A does not ex- 
cite the diffusive momentum density, which is generated 
by the transverse part of the source, but only excites the 
sound modes. The contribution to the medium excitation 
coming from the term proportional to A is important for 
generating Mach-like signals in the final azimuthal par- 
ticle spectrum. 

In what follows, we will use as our source term, 
treating A as an adjustable parameter and determining 
the rate of energy deposition dE/dt from different the- 
oretical models to be discussed in the next subsection. 
Our purpose here is not to suggest that the form of A 
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given in (f22j) is necessarily the correct form for the QGP 
created in heavy-ion collisions, but rather to motivate 
the general form of the source, (|2Tj) . The ansatz provides 
thus a connection between the hard, perturbativc QCD 
physics of jet quenching and the soft, nonperturbative 
QCD physics of medium response. 



D. Energy Deposition Scenarios 

In the following, we investigate three different scenar- 
ios for the energy deposition into the medium. In all 
three cases, we assume that the medium properties do 
not change as a function of space or time during the 
energy deposition. With this assumption, the spatiotcm- 
poral structure of the energy deposition dE/dt is inde- 
pendent of changes in medium properties and a function 
of the assumed physics of the source only. In addition, 
in order to study the dependence of the Shockwave signal 
on the functional form of dE/dt, we adjust the medium 
properties for each scenario such that the integrated en- 
ergy deposition AE to t = J^°dt dE/dt is the same. In 
contrast, if one would do calculations in a given micro- 
scopical model for parton-medium interactions and fix 
the medium to be e.g. at a given T, AEtot would not 
necessarily be the same in different scenarios of energy 
deposition. 

We adopt the procedure of normalizing to the same 
AEtot nevertheless because the relevant microscopical de- 
grees of freedom in the medium and their interaction with 
a hard probe are not known. In a more realistic model, 
with a hydrodynamically expanding medium, constraints 
from data on high Pt observables could be utilized in- 
stead by requiring each of the energy deposition scenar- 
ios to agree with the observed suppression of high Pt 
hadrons. In the absence of such data for a constant 
model, requiring the same AE to t is a substitute for such 
a constraint. 

In the first scenario, we assume that the source entering 
the medium is given by a single on-shell parton which 
interacts with the medium only elastically. The energy 
transfer into the medium is then given by the expression 

m 



Time Dependence of Energy Loss Rate 



a s C 2 m 2 D 2^/EpT 



m D 



(25) 



where a s — g 2 /(4ir) is the strong coupling, mrj = gT 
the Debye mass, C2 the appropriate color factor for a 
quark (4/3) or a gluon (3) and E p the energy of the hard 
parton. Under the assumption that the hard parton is 
sufficiently energetic such that E p 3> AE to t, the weak 
dependence on E p can be neglected and (^f ) c assumes 
a constant value if the medium temperature remains un- 
changed. Given this functional shape of dE/dt, we will 
refer to this scenario as Flat in the following. 

In a second scenario, we still assume that the source 
entering the medium is a single on-shell parton with suf- 
ficient energy such that E p 3> AE to t is realized, but we 
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FIG. 1: (Color online) The time dependence of the energy 
deposition rate into the medium for given AEtot = 20 GeV 
as calculated in three different scenarios Crescendo, Flat and 
Decreasing for the in-medium evolution of an initial hard par- 
ton (see text). 



also allow for inelastic interactions of this parton with 
the medium which induce radiation. As outlined in [f9j 
(see also [I(J), the gluons radiated from the hard par- 
ent partons become themselves sources of clastic energy 
loss with the medium and can also be absorbed by the 
medium if their energy becomes 0(T). 

In this formalism, the effect of radiated gluons to de- 
posit energy into the medium can be cast into the form 
of an evolution equation for the distribution function 
f(uj,t), which describes the distribution of radiated glu- 
ons in the medium at time t with energy uj. This evolu- 
tion equation reads 

P(^)-Il^)f(^ = ^ t (^) (^) 

where e(uj) is the collisional energy loss rate for gluons 
obtained from Eq. ([25]) as 



3 2 2VujT 

e(w) = -a s m D In 

2 mo 



(27) 



and in the Armesto-Salgado- Wiedemann (ASW) formal- 
ism [27I ] the spectrum of radiated gluons is given by 



dl 

dwdt 

where 



Vqa s C2 



[1 + i) tan 



Rc- 



,3/2 



q = 2a s C2'mf)T\rL 



2^/EpT 



(28) 



(29) 



is used to adjust the strength of the inelastic interactions. 
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Equation (|26|) must be solved numerically. Given such 
a solution, the rate of energy gained by the medium from 
radiated gluons is fl9l ] 



dl 



duje(uj) f (w, t) + / duiuj-, 
dt ) R L m , n Jo dudt (30) 

where Lj m in = T and u max = E p /2. The total energy 
deposition into the medium is then given as the sum of 
the collisional and radiative contributions 



dE _ (dE\ (dE 



dt 



dt J 



c 



dt 



(31) 



where (dE/dt)c is obtained from Eq. (|2"5|) . This second 
scenario leads to an increase of ^M- in time, therefore it 

at ' 

has been named Crescendo. 

In a third scenario, we take the source entering the 
medium to be a highly virtual parton which subsequently 
evolves into a parton shower. In addition, we do not make 
the assumption E p AE to t but consider finite energy 
kinematics for all partons. For this, we utilize the Monte 
Carlo (MC) code YaJEM (Yet another Jet Energy-loss 
Model) [HI, [2£| . In the following, we summarize the es- 
sential parts of the computation, details can be found in 
[29j (the scenario used in this paper corresponds to the 
DRAG (medium-induced drag force) scenario described 
in [11]). 

We model the evolution from the initial parton to a 
final state parton shower as a series of branching pro- 
cesses a — > b + c where a is called the parent parton 
and b and c are referred to as daughters. In QCD, the 
allowed branching processes are q — > qg, g — > gg and 
g — > qq. The kinematics of a branching is described in 
terms of the virtuality scale Q 2 and of the energy fraction 
z, where the energy of daughter b is given by Eb = zE a 
and of the daughter c by E c = (1 — z)E a . It is convenient 
to introduce t = lnQ 2 /Ag CjD where Aq C d is the scale 
parameter of QCD. t takes a role similar to a time in the 
evolution equations, as it describes the evolution from 
some high initial virtuality Qo (^o) to a lower virtuality 
Qm (t m ) at which the next branching occurs. In terms 
of the two variables, the differential probability dP a for 
a parton a to branch is [3(J Hl| 



dP a = ^P a ^bc(z)dtdz 



(32) 



where the splitting kernels P a ^bc{z) read 

(33) 
(34) 



1 + z 2 

P q ^ qg (z) = 



Pg^gg( z ) — 3 



{l-z{l-z)f 



P g ^(z) = N F /2(z 2 + (1 - zf). (35) 

We do not consider electromagnetic branchings. Np 
counts the number of active quark flavours for given vir- 
tuality. The resulting system of equations describing the 



branching processes in vacuum is solved numerically us- 
ing MC techniques utilizing the Pyshow code [32| . 

In order to make the link from momentum space where 
the shower evolution takes place to position space where 
the medium perturbations evolve, we assume that the 
average formation time of a shower parton with virtuality 
Q is developed on the timescale l/Q, i.e. the average 
lifetime of a virtual parton with virtuality Qb coming 
from a parent parton with virtuality Q a is in the rest 
frame of the original hard collision (the local rest frame 
of the medium may be different by a flow boost as the 
medium may not be static) given by 



Eb Eb 



(36) 



We assume that the actual formation time can then be 
obtained from a probability distribution 



P(Tb) = cxp 



n 
(n) 



(37) 



which we sample to determine the actual formation time 
of the fluctuation in each branching. 

We assume that the medium induces an approximately 
constant energy loss on each propagating parton. The 
medium is then characterized by a drag coefficient D 
which describes the energy loss per unit pathlcngth. In 
the simulation, the energy (and momentum) of each 
propagating parton are reduced by 



d(D 



(38) 



For a propagating gluon the energy loss is increased by 
the color factor ratio 2.25. While Eq. (|3"5)) describes 
the mean energy loss, the actual energy loss due to the 
medium is randomized in each event. 

The dynamics of energy deposition in YaJEM is rather 
different from the Crescendo scenario. The initial branch- 
ing processes down from a highly virtual state happen 
very fast and lead to a pronounced initial rise in energy 
deposition as the number of partons undergoing elastic 
energy loss increases. However, the finite energy of the 
parton shower which is explicitly considered leads to a 
turnover: As partons become absorbed by the medium, 
the number of additionally radiated gluons is limited by 
kinematic constraints. As a result, the functional shape 
of 4f is decreasing in time, therefore the scenario is la- 
beled Decreasing in the following. 



III. RESULTS 

A. Input Values and Qualitative Expectations 

In the results below we consider the following situ- 
ation: a source parton is created at time t = and 
travels through the medium until the energy deposition 
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ceases at the time t = 5 fm/c. As discussed in section 
III C[ the source parton is assumed to excite the medium 
through the source term ([21]) . We continue to evolve the 
medium response until time t = 7 fm at which point the 
medium is assumed to hadronize into the spectrum given 
by (fT2| . Performing the freeze-out 2 fm after the source 
is turned off improves the validity of the linearized hy- 
drodynamic assumption because the peak amplitude of 
the energy and momentum density perturbations decays 
in time once the source is turned off. 

The underlying medium is at temperature To = 250 
MeV and has speed of sound c 2 = 1/3. The shear viscos- 
ity to entropy density ratio, rj/s, is treated as an input 
parameter that we vary from r]/s = 0.1 — 0.2. The range 
of input values for r\js are consistent with phenomeno- 
logical observations from heavy-ion collisions [33j | . The 
local excitation parameter, A, mentioned in section III CI 
is treated as an input parameter, which we vary between 
— 2 fm. Although we treat A and 77/s as parameters, 
rigorous determination of their values must come from 
the underlying theory. One may ask what sets an upper 
limit for A. From the point of view of hydrodynamics, A 
appears as the coefficient of a gradient, so in principle its 
value sets a minimum resolution for the hydrodynamic 
equations of motion. However, in the case we are consid- 
ering, we allow the medium response to decay in time, 
improving the validity of hydrodynamics, and in particu- 
lar, linearized hydrodynamics, once the freeze-out occurs. 
To quantitatively set an upper limit on A, one must solve 
for the medium response for a given scenario and see how 
large the perturbations are. In each of the scenarios we 
consider below, we find that the linearized approximation 
holds (to the extent that 8u appearing in (|lip remains 
less than 1). 

The primary goal of this paper is to demonstrate 
that the time dependence of the energy deposition rate, 
dE/dt, which appears as a coefficient to the source (|2"Tj) 
is highly relevant for the appearance of a Shockwave sig- 
nal in the azimuthal spectrum (|12p . We will consider 
three different energy deposition scenarios, Crescendo, 
Flat and Decreasing, that have been discussed in detail 
in section Hi Dl In order to make the comparison consis- 
tent, in all three scenarios the total energy deposited into 
the medium is 20 GeV (see Figure [1]). Before examining 
the results of our calculations it is useful to qualitatively 
consider what may be expected. From equations (|lip 
and (p~2|) one can see that the appearance of a Shockwave 
signal is sensitive to the combination 5u{x) pt/T, where, 
as discussed in section Hi Bl Su(x) is the four- velocity in- 
duced by the source. In the limit of validity of linearized 
hydrodynamics, Su < 1, thus, as mentioned in a previous 
work |34f . one expects the Shockwave signal to be more 
pronounced for larger values of pt/T. However, it's not 
hard to see that this expectation remains even for non- 
linear hydrodynamics. Even if one doesn't linearize the 
disturbance created by the hard parton, equations (TTTj) 
and (|12p still predict an enhanced signal for larger px- 

Experimentally, the opposite trend is seen 0, Q : the 



double-peaked structure observed in the away-side distri- 
bution of di-hadron correlations is more pronounced for 
smaller values of px- This observation does not mean the 
Mach cone Shockwave is not responsible for the double- 
peaked structure seen experimentally, but rather exposes 
the limitations of hydrodynamics. As px increases, the 
driving mechanism behind the correlation structure shifts 
from bulk recoil to hard fragmentation. Most events pro- 
duce correlations in the low px hydrodynamical regime, 
but at higher px, there's an increased bias to see the 
comparatively rare events in which hard pQCD is the 
mechanism underlying the correlation. At px beyond 
5-6 GeV, these hard events completely dominate the vis- 
ible correlation. It is this transition from soft to hard 
physics which governs the transition from shoulder-region 
to head-region in the data and is beyond the scope of our 
current work (see [35| for an analysis of this transition). 

Generally speaking, one would expect that larger am- 
plitudes of the induced energy and momentum density 
perturbations would be more likely to generate a shock- 
wave signal in the azimuthal spectrum than smaller ones. 
In terms of the time dependence of the energy deposition 
rate, dE/dt, the naive expectation is that an energy de- 
position rate that grows in time is more favorable for 
the appearance of a Shockwave signal than one which de- 
creases or is flat. To make the argument more concrete, 
consider a source that deposits some amount of energy 
and momentum at an instant in time and is then turned 
off. The time dependence of the medium perturbation 
due to such a source can be tracked by going back to ([3]- 
[5]) and using the reverse Fourier transform © . Perform- 
ing the uj integration using contour methods and leaving 
the k integration undone shows that, to first order in 
shear viscosity, as a function of time 



Se, g L ~ exp 



and 



gx ~ exp 



3T,k' 2 t 



(39) 



(40) 



The medium excitation decays exponentially as a func- 
tion of time, meaning that the azimuthal spectrum (|12[) 
more reflects the strength of energy deposition at later 
times than earlier times. However, the situation is com- 
plicated by the fact that the diffusive contribution, gx, 
decays more quickly in time than the sound contribu- 
tion, Se and gL- It was mentioned in section [II CI that 
the appearance of a Shockwave signal in the azimuthal 
spectrum will come from 5e and gL , whereas gx tends to 
fill up any double-peak structure. Even though the am- 
plitude of the medium excitation decays exponentially in 
time for both the sound and diffusion modes, the ratio of 
the sound/diffusive contribution actually grows in time. 
From these considerations, it is not immediately clear 
whether an energy deposition rate which grows in time 
or decreases in time is more favorable for generating a 



8 




FIG. 2: (Color online) The azimuthal hadron spectrum (I12|l for the case of n/s = 0.10 and local excitation parameter, A = 
fm. The shapes of the spectrums for the different energy deposition scenarios are essentially the same for these parameters. 
The magnitude of the signal is larger for Crescendo , which grows in time, than for Flat or Decreasing , even though the same 
total energy is deposited in each case. The conical structure generated by the source (|2ip when A = fm is not strong enough 
to overcome the diffusive wake in the final spectrums, where one sees a single peak at cj> — 0, which defines the direction of 
source propagation. The larger (blue) tick marks on the cj> axis indicate where one would naively expect conical peaks to appear 
for the speed of sound used here, <j> — arccos c s . 




FIG. 3: (Color online) The same as in Figure [2] but now A = 0.5 fm. The spectrum is noticeably flatter than seen in the A = 
fm case, especially for the pr — 3 — 4 GeV plot, where a double-peak begins to emerge in the Crescendo curve. The difference 
in the shapes of the spectrums generated by the different energy deposition scenarios is most noticeable when comparing the 
Crescendo and Decreasing scenarios. The energy deposition which grows in time generates a stronger conical signal. The larger 
(blue) tick marks on the <j> ax i s indicate where one would naively expect conical peaks to appear for the speed of sound we have 
used. 



shockwavc signal. In the next subsection we present the 
results of our calculations. 



B. Numerical Results 

In this section we present results for the azimuthal 
hadron spectrum (|12|) obtained using the medium param- 



eters and energy deposition scenarios discussed above. In 
all results we will show the spectrum for the three en- 
ergy deposition scenarios, Crescendo, Flat and Decreas- 
ing, and for the bins px = 2 — 3 and pr = 3 — 4 GeV. We 
now present results for the azimuthal hadron spectrum 
(|12[) obtained using the medium parameters and energy 
deposition scenarios discussed above. In all results we 
will show the spectrum for the three energy deposition 
scenarios, Crescendo, Flat and Decreasing, and for the 



9 




FIG. 4: (Color online) The same as in Figures [2] and [3] but now A = 1.0 fm. The double-peaked structure is significantly 
enhanced, particularly in the Crescendo and Flat spectrums. The difference in the shapes of the spectrums generated by the 
different energy deposition scenarios is obvious. The larger (blue) tick marks on the <j> axis indicate where one would naively 
expect conical peaks to appear for the speed of sound used here. 




FIG. 5: (Color online) The azimuthal hadron spectrum (|12[) for the case of n/s = 0.15 and A = 1.0 fm. The effect of viscosity to 
smear out any double-peaked structure is obvious when comparing to the results in Figure U where we considered n/s = 0.10 
and A = 1.0 fm. The larger (blue) tick marks on the 4> axis indicate where one would naively expect conical peaks to appear 
for the speed of sound used here. 



bins pt = 2 — 3 and pr = 3 — 4 GeV. Results for shear 
viscosity to entropy density ratio rj/s = 0.10 and local 
excitation parameter, A = fm (A is discussed in section 
III C[) are shown in Figure [2] Recall that the direction of 
the source propagation determines the direction of <fi = 
in our plots. The results show that the shape of the 
spectrum is roughly the same for all three energy depo- 
sition scenarios, however, the magnitude of the signal is 
larger for Crescendo, which grows in time, than for Flat 
or Decreasing. This change in magnitude between the 
different energy deposition scenarios reflects the viscous 
nature of the medium. Energy which is deposited at ear- 



lier times (such as in the Decreasing scenario) has more 
time to equilibrate with the background medium before 
freeze-out. What is noticeably missing in Figure [5] is the 
appearance of a double-peaked structure, or Shockwave 
signal. Apparently the conical structure generated by 
the source (f2"Tj) when A = fm is not strong enough to 
overcome the diffusive wake in the final spectrum. 

Next we consider results for rj/s = 0.10 and A = 0.5 
fm which are shown in Figure [3] The spectrum is notice- 
ably flatter than seen in the A = f m case, especially for 
the pt — 3 — 4 GeV plot, where a double-peak begins to 
emerge in the Crescendo curve. Notice that not only are 
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FIG. 6: (Color online) The azimuthal hadron spectrum for n/s = 0.15 and A = 1.5 fm. Increasing A restores the Shockwave 
signal in the Flat and Crescendo energy deposition scenarios, (compare to Figure [5] where A = 1.0). However, consistent with 
the above results, the Decreasing energy deposition remains mostly flat. The larger (blue) tick marks on the (j> axis indicate 
where one would naively expect conical peaks to appear for the speed of sound used here. 




FIG. 7: (Color online) The resulting spectrums for n/s — 0.20 and A = 2.0 fm. Again, the result of increasing the viscosity 
is to smear out the double-peaked structure, even though we have increased A by the same fraction as presented in Figure [5] 
Recall that the larger (blue) tick marks on the <j> axis indicate where one would naively expect conical peaks to appear for the 
speed of sound used here. 



the magnitudes of the curves for the different energy de- 
position scenarios different, but also the shapes. The dif- 
ference is most noticeable when comparing the Crescendo 
and Decreasing scenarios. The shape of the Decreasing 
spectrum is very similar to the case of A = fm, how- 
ever the shape of the Crescendo spectrum is much flatter, 
even showing a slight dip at <j) = 0. The energy deposi- 
tion which grows in time appears to generate a cleaner 
signal of the underlying conical structure induced in the 
medium. 

We continue by considering results for 77/s = 0.10 and 
A = 1.0 fm presented in Figure 21 Here the double- 



peaked structure becomes quite pronounced for the Flat 
and Crescendo curves, however the Decreasing result re- 
mains mostly flat, with a slight dip at <fi = in the 
Pt = 3 — 4 GcV range. The Flat and Crescendo spec- 
trums have a similar shape, however the Crescendo result 
shows a moderately more pronounced double-peak. The 
trend is emerging that for fixed values of 77/s and in- 
creasing values of A, an energy deposition scenario which 
increases in time generates a more pronounced Shockwave 
signal in the final azimuthal spectrum. 

The results for 77/s = 0.15 and A = 1.0 fm are shown 
in Figure [5] One can see that a 50 percent increase in 
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P T = 2-3 [GeV], X = 1 .25 [fin] ( a ) P T = 2-3 [GeV], A. = 1 .25 [fin] (b) 

n/s = 0.10 rj/s = 0.10 




4> <l> 

FIG. 8: (Color online) Results for pt = 2 — 3 GeV, 77/s = 0.10 and A = 1.25 fm. The left panel shows the result without 
including the boundary and derivative terms (as discussed in the text), whereas the right panel contains those contributions. 
In the right panel, one sees that the difference in the shapes of the curves are less significant, and also narrower. This results 
from inward flow created by the absorption of the source at an instant in time, and is an artifact of an infinite energy (for the 
case of the Crescendo and Flat scenarios), static medium assumption. We have not included the boundary terms in the above 
results, as they obscure the physics we are trying to study. The larger (blue) tick marks on the <j> axis indicate where one would 
naively expect conical peaks to appear for the speed of sound used here. 



the shear viscosity has a significant effect on the final 
azimuthal spectrum by comparing with Figure UJ The 
double-peaked Shockwave signature has been smeared 
out by the viscous effects. Continuing with rj/s = 0.15 
and increasing A to 1.5 fm restores the Shockwave signal 
in the Flat and Crescendo energy deposition scenarios, as 
is evident from Figure [6] However, consistent with the 
above results, the Decreasing energy deposition remains 
mostly flat. Finally, we present results for rj/s = 0.20 and 
A = 2.0 fm in Figure [71 Again, the result of increasing 
the viscosity is to smear out the double-peaked structure, 
even though we have increased A by the same fraction. 

We here discuss the form of the source, which 
we have employed for the results presented above. As 
mentioned briefly in section II C, we do not include a 
time derivative on dE/dt or boundary terms which arise 
from the source being turned on and off at t = and 
t = T. As noted above, we have explicitly checked that 
the derivative term is numerically insignificant for the 
energy deposition scenarios we consider and can safely be 
ignored. However, the boundary terms, in particular, the 
term resulting from turning the source off at time t = T, 
is not numerically insignificant. In fact, the effect of the 
boundary term at t = T is to create a strong inward flow 
as the source is absorbed by the medium. This inward 
flow tends to destroy the conical Mach cone signal, and 
reduce the differences in the three scenarios. 

This feature is demonstrated in Figure [8] where we 
show the azimuthal spectrum for the bin px = 2 — 3 
GeV, rj/s = 0.10 and A = 1.25 fm. The left panel shows 
the result without including the boundary and derivative 



terms, whereas the right panel contains those contribu- 
tions (these are dominated by the t = T boundary piece). 
The left panel is characteristic of the plots shown above, 
in that the Crescendo curve has the most pronounced 
double peak structure, whereas the Flat and Decreasing 
curves are less pronounced. In the right panel, however, 
one sees that the difference in the shapes of the curves 
are less significant, and also narrower. This results di- 
rectly from inward flow created by the absorption of the 
source at an instant in time. 

The absorption at t — T is an artifact of our proce- 
dure to stop the simulation, and tends to obscure the 
physics we are trying to study. In a real physical sit- 
uation, the source is not simply absorbed at some fixed 
instant, but may escape the medium or run out of energy, 
etc. With an infinite energy (for the case of the Crescendo 
and Flat scenarios), static medium assumption, the prob- 
lem of stopping the energy deposition is not well defined 
(it would go on forever), so we have to make a choice 
as to what we want to show as a freeze-out distribution. 
We believe we can get closer to what we're interested in 
by dropping the boundary terms - which we have done 
in the above. However, the complications just described 
illustrate the need for a more realistic description of both 
the medium evolution and the source term. 

We conclude this section by briefly considering how our 
results depend on the specific choice of speed of sound 
c 2 s = 1/3 and temperature T = 250 MeV we have used 
here. Both the speed of sound and the temperature ap- 
pear explicitly in the medium distribution function (sec 
© and (TTTj) ) used for freeze-out. Changing the temper- 
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P T = 2-3 [GeV], X = 1 [fm] 
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FIG. 9: (Color online) Results for p T = 2-3 GeV, v/s = 0.10 
and A = 1.0 fm with two different speeds of sound. As dis- 
cussed in the text, one would anticipate the spectrum be- 
comes more pronounced when lowering the speed of sound 
and that the peaks would appear at larger opening angles. 
These two features are indeed observed. However, notice that 
the peak angles do not correspond to a naive prediction based 
on purely geometric arguments (the blue ticks in the Figure 
would naively correspond to c s = 0.577) even for the simple 
scenario considered here. 



ature in the medium distribution function (|llj) will have 
an effect on the strength of the freeze-out signal. Specif- 
ically, lowering the temperature will tend to generate a 
stronger signal. Likewise, lowering the speed of sound in 
the medium distribution function will also tend to create 
a stronger freeze-out signal (even more so than the tem- 
perature, since the speed of sound appears only with the 
momentum flow). 

The temperature also enters the sound attenuation 
length, r s = g^r, which appears in the equations of mo- 
tion for linearized hydrodynamics (J3] - [5J) - Changing the 
temperature in T s is effectively like changing rj/s, which 
has been analyzed in the results above, thus we do not 
consider that aspect further here. However, the speed of 
sound appears in a non-trivial way in the equations of 
motion. As is well understood, the speed of sound gov- 
erns the angle of propagation of Mach cone shock waves 
generated by a projectile. We present in Figure [9] the ef- 
fects of changing the speed of sound on the the azimuthal 
hadron spectrum for the Crescendo and Flat energy de- 
position scenarios and for the bin pt = 2 — 3 GeV. The 
results arc for rj/s = 0.10, A = 1 fm and speed of sound 
c s = 0.577 and c s = 0.5. As anticipated, the spectrum 
becomes more pronounced when lowering the speed of 
sound, and the peaks also appear at larger opening an- 
gles. However, notice that the peak angles do not cor- 
respond to a naive prediction based on purely geometric 
arguments even for the simple scenario considered here 
(the blue ticks in the Figure would naively correspond 



to c s — 0.577). One thus needs to be very careful not 
to interpret the experimentally measured opening angle 
of the correlation geometrically as directly related to the 
speed of sound — in a realistic medium, the combination 
of trigger bias with a longitudinal, transverse and ellip- 
tic flow field will have an even more significant influence 
on the angular structure than in our simplified medium 
study. 

IV. CONCLUSIONS 

We have conducted a systematic study of the shock- 
wave excitation for different models of the spatio- 
temporal structure of energy deposition into the medium 
within a linearized hydrodynamical framework in a con- 
stant medium under different assumptions with regard to 
medium properties and the interaction of the source with 
the medium. The results exhibit a few generic trends: 

• Only a strong gradient term (A > 0) in the source 
(|2"T|) leads to an observable double-hump structure. 
The observation that a gradient term is necessary 
to excite an observable double-peak has been made 
early on (Toj . however, our implementation of the 
source (|21[) provides a way to quantify how strong 
the gradient term must be. 

• For fixed strength of the gradient term, viscous ef- 
fects (i.e. larger values of rj/s) weaken the double- 
hump structure. This has a natural explanation in 
terms of entropy generation dissipating the shock- 
wave, but even for relatively small rj/s the effect 
appears rather pronounced. 

• Consistently for all assumptions about the struc- 
ture of the source term and the medium shear vis- 
cosity, an energy deposition dE/dt which increases 
as a function of time leads to more pronounced 
shockwavc-like correlations than dE/dt decreasing 
in time — both in the absolute strength of the cor- 
relation as well as in the shape. This effect is not 
small — the correlation strength can be reduced 
more than an order of magnitude in the 'Decreas- 
ing' as compared to the 'Crescendo' scenario, al- 
though the precise factor depends on the medium 
properties. 

Applied to the Shockwave interpretation of measured 
correlations in heavy-ion collisions, these findings imply 
that the measured signal strongly depend on medium 
properties, the local structure of the source term and the 
spatio-temporal pattern of energy deposition. Thus, if 
one is able to model the evolution of an initial hard par- 
ton and its interaction with the medium with sufficient 
precision, one can determine medium parameters like rj/s 
from the measured correlations. On the other hand, if 
one can extract the medium properties with other meth- 
ods, the correlations then place tight constraints on the 
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dynamics of hard parton evolution in the medium. Which 
road will be taken first still remains to be seen. 
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